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We numerically investigate statistical ensembles for the occupations of eigenstates of an isolated 
quantum system emerging as a result of quantum quenches. The systems investigated are sparse 
random matrix Hamiltonians and disordered lattices. In the former case, the quench consists of 
sudden switching-on the off-diagonal elements of the Hamiltonian. In the latter case, it is sud- 
den switching-on of the hopping between adjacent lattice sites. The quench-induced ensembles are 
compared with the so-called "quantum micro-canonical" (QMC) ensemble describing quantum su- 
perpositions with fixed energy expectation values. Our main finding is that quantum quenches with 
sparse random matrices having one special diagonal element lead to the condensation phenomenon 
predicted for the QMC ensemble. Away from the QMC condensation regime, the overall agreement 
with the QMC predictions is only qualitative for both random matrices and disordered lattices but 
with some cases of a very good quantitative agreement. In the case of disordered lattices, the QMC 
ensemble can be used to estimate the probability of finding a particle in a localized or delocalized 
eigenstate. 

PACS numbers: 



INTRODUCTION 



The issue of statistical ensembles emerging in isolated 
quantum systems after a sudden perturbation, often re- 
ferred to as "quantum quench" , has recently come into 
a sharp focus of statistical physics research Of par- 
ticular interest for the present paper are the statistical 
ensembles for the occupations of the energy eigenstates 
of an isolated quantum system. We consider two kinds 
of quenches: a more abstract one induced by switching- 
on the off-diagonal elements of a sparse random ma- 
trix Hamiltonian, and a more physical one induced by 
switching-on random hopping elements in a disordered 
potential. Our interest was motivated by the numeri- 
cal result reported in Ref. that, for a certain class 
of random matrices, the above quench induces the so- 
called "quantum micro-canonical" (QMC) ensemble 0- 
[T3 |. The QMC ensemble, when treated analytically, 
exhibits a remarkable phenomenon of condensation ob- 
tained in Rcfs. [l]l G3- However, the numerical in- 
vestigation of Ref. Q could not reach the condensation 
regime with the class of random matrices considered. 

In this paper, we investigate a new class of random ma- 
trices with a special diagonal element, which allows us to 
reach the QMC condensation regime. At the same time, 
we show that the emerging ensembles often show notice- 



able differences from the QMC predictions. We find that 
the QMC condensation is connected to the localization 
properties of the eigenstates near the edge of the energy 
spectrum. We also find that the ensembles emerging in 
quenched disordered potentials are less universal than in 
the case of quenched random matrices. However, they 
are still in a qualitative agreement with the QMC predic- 
tions. The random matrix quench may possibly be used 
to mimic the effect of multiple small quenches on phys- 
ical systems with a large number of quantum levels but 
not too large a number of particles |6j. The disordered- 
potential quench is relevant to the recent experiments (l5r- 
llal . which report the evidence of Anderson localization 
[20(. It can, in particular, help to estimate the probabil- 
ity of a particle being below or above the mobility edge 
after the quench. 

The rest of the paper is organized as follows: In Sec.UTl 
we give the general formulation of the problem. In Sec. 
IIII1 we describe the QMC ensemble. Sections IIVI and 
IVl deal with the random matrix Hamiltonians and the 
disordered potentials, respectively. 



II. GENERAL FORMULATION 

We consider an isolated quantum system having a large 
but finite number of quantum states N with a Hamilto- 
nian that can be decomposed as follows: 
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In the context of sparse random matrices, % represents 
the entire random matrix, Hq the diagonal elements of 
H and Hpcrt the off-diagonal elements. In the case of 
disordered lattices, T-Lq represents the on-site energies, 
and Hpcrt the hopping terms. We denote the eigenstates 
of % as ipi and the corresponding eigenenergies as Ei. 
The eigenstates of the Hamiltonian Hq are to be denoted 
as (f>k . We also choose the origin of the energy axis such 
that Yli=i Ei = 0, i.e., for a symmetric energy spectrum, 
Ej = corresponds to the middle of the spectrum. 

We ask the following question: if the system is initially 
in an eigenstate <pk of the Hamiltonian and then 
"Hport is switched on, what kind of statistical ensemble 
emerges for the occupations of the eigenstates ipi of the 
Hamiltonian T~Cl Specifically, we consider the expansion 



(2) 



where Ci are the complex amplitudes, and investigate the 
average values for the occupations of eigenstates (|ci| 2 ) 
as a function of energy Ei and energy E av = {4>k\H\4>k) = 
(4>k\'Ho\4>k) ■ The latter equality is due to the fact, that, in 
the both cases considered, % pe rt does not have diagonal 
elements in the eigenbasis of Ho- In other words, the 
fc-th diagonal element of Ho fixes the energy expectation 
value of the of <pk with respect to H. For convenience, we 
also introduce the "occupations of eigenstates" pi = \ci\ 2 . 
The averaged occupations of eigenstates as a function of 
their energy is to be denoted as (pi)(Ei). In the following, 
we often drop index i in this function and represent it as 
(P)(E). 

The quantity (pi)(Ei) is to be obtained numerically by 
organizing eigenstates ipi with close values of Ei into "en- 
ergy bins" and then by averaging within each bin. The 
statistics will be further improved by including within a 
given energy bin all piS for the ensemble of initial states 
4>k with close values of E av and for different random re- 
alizations of random matrix Hamiltonians or disordered 
potentials. 

The present work is complementary to a number of 
studies that considered the statistics of eigenvectors for 
band random matrices (2l| - [23| and quantum quenches for 
other special kinds of sparse random matrices (8[. The 
distinctive focus here is on the statistics for a fixed value 
of E av . If one were to perform the above averaging over 
all initial states i.e. not just over those having close 
values of E av , this would result in (p)(E) = 1/N, i.e. 
there would be no energy dependence of (p). However, 
the constraint of the energy expectation value E av is in- 
compatible with the constant value of (p). For example, 
if E av is below the middle of the energy spectrum of "H, 
then it should be expected that (p)(E) is a decreasing 
function of E. In general, this decrease is not expected 
to be exponential. Investigation of the character of this 
decrease is the goal of the present paper. 



III. THE QMC ENSEMBLE 

The QMC ensemble is a quantum ensemble for which 
the constraint on the energy expectation value is the only 
constraint in the Hilbert space of the system d, [ll|, [l3| . 
In the Euclidean space of variables pi, the QMC ensem- 
ble is defined by the uniform probability density on the 
manifold constrained by the conditions 
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EiPi = E a 
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Pi = 1, 



and 



Pi > 0. 



(3) 



(4) 



(5) 



The phases of the complex amplitudes Cj are completely 
random. 

Contrary to naive expectations, the QMC ensemble 
does not lead to the conventional thermodynamics. The 
marginal probability distribution Pi(pi) associated with 
the QMC ensemble was derived in Ref. @. In the most 
common case of pi <SC 1 , the result is 



Piipi) = Pl (0)e- N ^ 1+X ^- E ^\ 



(P, 



N[l + \(E t -E av )Y 
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The parameter A has the meaning of the inverse Hilbert 
space temperature. It can be calculated numerically by 
substituting the averaged probabilities ([7]) into the con- 
straint © on the energy expectation value: 



N 



0. 



(8) 



In Refs. 0, [TH it was also shown that, for macroscopic 
systems consisting of N s weakly interacting constituents 
and for E av < - 0(Ei/y/{N a ), the QMC ensemble 
implies condensation into the ground state, which means 
that the ground state acquires large average occupation 



(Pi) « ^ 



(9) 



(typically of the order of 1), while all other eigenstates 
have occupations of the order of 1/N given by formula 
([7]). The probability distribution P\(pi) in this case is 
not given by formula ©• Instead, it is narrowly peaked 
around pi = (pi). 

For other kinds of energy spectra or for non- 
macroscopic systems, the QMC condensation is also pos- 
sible, but it may be accompanied by the anomalously 
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large occupations of a few lowest energy levels instead of 
just one level[ll[I|. 

The above condensation together with the very slowly 
decaying tail of form (O leads to a non-thermal den- 
sity matrix for a small subsystem within a macroscopic 
system[|[3[Il]. 

In the following, the occupation numbers (pi) gener- 
ated as a result of numerically simulated quenches are 
compared with prediction (O for the QMC ensemble. In 
order to obtain A, we combine Eq. (|7|) and Eq. ©. The 
result is 

y Ei ~ Eav = o. (io) 

The value of A is then obtained numerically by locating 
the zero of the LHS of Eq. (TU]) closest to A = 0. 

Below, we also compare the quench-generated ensem- 
bles with the canonical thermal ensemble corresponding 
to temperature T, which is chosen such that the thermal 
average energy is equal to E av : 

1 N 

E av = -Y,Eie-^. (11) 

i=l 

Here Z is the partition function. In principle, there is no 
fundamental reason for the thermal ensemble to emerge, 
but such a comparison still provides a useful perspective. 

For the random-matrix quench considered in this work, 
there is an argument why the QMC ensemble may be 
expected to emerge. Namely, the transformation between 
the eigenbasis of the Hamiltonian "Ho and the eigenbasis 
of the Hamiltonian H can be reasonably expected to be 
close to a random rotation in the Hilbcrt space. If it 
were indeed a random rotation, and then, in addition, the 
constraint on the energy expectation value were imposed, 
then the resulting set of states <frk would indeed represent 
the QMC ensemble. 

In general, however, for physically realizable quenches, 
the Hamiltonian % has the character of a band matrix in 
the basis of the Hamiltonian Hq. Therefore, the trans- 
formation between the eigenbasis of H and the eigenbasis 
of Hq is not supposed to have the character of a random 
rotation. Yet, if the system is not too large, then one can 
expect that the QMC ensemble would provide at least a 
rough guide to the behavior of the function (p) (E) . In 
particular, it was shown in Ref. @ that multiple small 
quenches applied to finite clusters of spins 1/2 lead to 
QMC-like ensembles. If a random matrix quench also 
produces a QMC-like ensemble, it may imply that the 
effect of multiple small physical quenches may be repre- 
sented by a single random-matrix quench. 

IV. RANDOM MATRIX HAMILTONIANS 

All numerical results to be presented in this section are 
obtained for random matrices of size NxN = 4096 x 4096 




FIG. 1: Normalized densities of states v(E) for two 4096 x 
4096 sparse random matrix Hamiltonians H described in the 
text: thick black line — all non-zero elements are real number 
randomly selected from the interval [—1,1]; thin grey line — 
the same Hamiltonian but with one special diagonal element 
equal to -4. 



with all diagonal elements and a fraction 30/4096 of the 
off-diagonal elements being different from zero. The non- 
zero entries are to be selected from different ensembles. 

The random matrix quench for the matrix of the above 
size with all non-zero entries being real and chosen ran- 
domly from interval [— 1 , 1] was first investigated (with- 
out calling it a "quench") in Ref. The density of 
states for this random matrix is shown in Fig. [T] The 
resulting occupations of eigenstates appeared to be very 
well describable by the QMC formula (JT)). However, as 
we illustrate below, the very good agreement reported 
in Ref. Q was somewhat accidental. As we discovered 
later, the statistical fluctuations between different initial 
states and different random matrix realizations are larger 
than the visible statistical noise for the occupations of 
energy bins. More complete averaging confirms overall 
satisfactory agreement, but also reveals noticeable dis- 
crepancies between the numerically generated ensembles 
and the QMC prediction, in particular, near the edges of 
the energy spectrum. For example, Fig. [3J shows these 
discrepancies for one of the values of E av used in Ref. . 
We return to the implications of Fig. [5] later in subsec- 
tion hvb 

Here we focus on the fact that the interesting regime 
of possible QMC condensation could not be reached with 
the above random matrix ensemble because of the follow- 
ing constraint: The QMC condensation regime requires 
E av to be sufficiently close to the the ground state energy. 
In this case, however, the lowest possible energy expec- 
tation value E av for the initial state <f>k is the smallest 
possible value of the diagonal element of H, i.e. — 1, 
which is not close enough to the ground state energy 
equal to about —6 [see Fig. Q]. Making the random 
matrix larger or less sparse would only lower the ground 
state energy without changing the minimum value of E av . 
On the other hand, making the matrix more sparse can 
indeed bring the ground state energy closer to the mini- 
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FIG. 2: (Color online) Occupations of eigenstates after a ran- 
dom matrix quench with all non-zero matrix elements in H be- 
ing real and selected randomly from the interval [—1,1]. The 
energy expectation value for the initial state is E av = —0.9. 
Numerical results represent 100 independent realizations of 
random matrices. The occupations of eigenstates (p) are aver- 
aged within energy bins, each containing N\, = 100-32 = 3200 
states. The bin averages are computed according to formulas 
lfl2) and . The QMC prediction is represented by red 
solid line. The exponential thermal distribution with same 
average energy is represented by the dashed green line. 



element as the initial state before the quench, i.e. E av 
becomes smaller than -1. 

We start with the example of E av = —2 illustrated 
in detail in Fig. 03 Figure 02 a ) shows unaveraged oc- 
cupations of eigenstates after the quench for a single 
realization of the random matrix Hamiltonian %. Fig- 
ures G3Jb) and (c) show the occupations of eigenstates 
averaged over 20 realizations of random matrices. All 
pairs {Ei,pi} from all 20 realizations are combined in 
an energy-ordered list, and then divided into the energy 
bins, with each bin containing N\, entries. Fixing the 
same N\, for each bin implies that the bins acquire vary- 
ing energy widths and may become particularly broad 
near the edges of the energy spectrum. Here and ev- 
erywhere below, the average characteristics for the j-th. 
energy bin are computed as follows: 



E, 



E Pi E i 
j-th bin 

E Pi 

j-th. bin 



Pi 



(12) 
(13) 



j-th bin 



mum value of E av , but, at the same time, the basis states 
4>k become too correlated with the eigenstates ipi . In the 
limit, when all off-diagonal elements are equal to zero, 
the eigenstates of Hq trivially coincide with the eigen- 
states of H. Therefore, on the approach to this limit, the 
random rotation argument given in Section Mil becomes 
increasingly inadequate, and the ensemble generated by 
the quench becomes very different from QMC. 

The main part of the next subsection introduces and 
investigates the new ensemble of random matrices con- 
taining one special diagonal element selected differently 
than the rest. This allows us to realize the QMC conden- 
sation regime. The idea is that, on the one hand, only one 
special diagonal element does not strongly influence the 
shape of the energy spectrum, while, on the other hand, 
this element gives us a predetermined value of E av which 
can be chosen much closer to the ground-state energy of 
rl . We further investigate ensembles with more than one 
special diagonal element. Finally, we will show that the 
result of Ref. Q extends to a complex- valued random 
matrix selected from the Gaussian unitary ensemble. 

A. Results 

1. Ensembles with one or more special diagonal elements 

In this subsection, we consider sparse random matri- 
ces, where all non-zero off-diagonal elements and all but 
one diagonal elements are chosen randomly from a uni- 
form distribution in the interval [-1,1]. A single diagonal 
element is assigned a value smaller than — 1. We then 
choose the basis state 4>k corresponding to this special 



The difference between Figs. G3Jb) and (c) is that, in the 
former case, the bin size is smaller, N\, = 20 • 32 = 640, 
while in the later case it is larger, TVb = 20 • 128 = 2560. 
Finally Fig. EJd) shows the inverse participation ratio 
(IPR) for the expansion of eigenstates ip, in terms of the 
basis states 4>t — see Eq. (fT4"|) . 

As seen in Fig. [3jb) the bin-averaged occupations for 
all but the lowest energy bin follow the QMC statistics 
very well. It is also very clear that the occupations of the 
eigenstates do not decrease exponentially as one would 
expect from a thermal distribution but rather decrease 
according to the power law of the QMC ensemble. In 
Fig. 03c), the discrepancy for the lowest energy bin was 
eliminated by choosing the bin size such that the lowest 
energy bin includes all the eigenstates with anomalously 
high inverse participation ratio. This aspect is to be dis- 
cussed further in subsection IIVBI 

Figure 03 a ) indicates that a quantum quench with 
i? av = — 2 does not exhibit the condensation into a single 
cigenstate, yet. However, by reducing the value of the 
special diagonal clement further, the QMC condensation 
can indeed be reached. In Fig. 0] we present the result 
of the quench for the same kind of random matrix but 
with the special diagonal element equal to -4. Fig. Hta) 
shows that the ground state has occupation comparable 
to 1, while the occupations of all other states are very 
small and, as shown in Fig. SJb), accurately describable 
by the QMC statistics. As shown in Fig. [1] the density 
of states in this case is still nearly the same as in the case 
of random matrix without a special diagonal element. 

In Fig. 02 we plot the occupation p\ of the ground state 
as a function of E av . The numerically generated values 
of pi strongly fluctuate as a function of i? av , because 
each plot point was generated with different realization 
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FIG. 3: (Color online) Occupations of eigenstates after a ran- 
dom matrix quench for the sparse Hamiltonian H with one 
special diagonal element as described in the text. The en- 
ergy expectation value for the initial state is E av = —2. (a) 
Unaveraged occupations of eigenstates p for a single realiza- 
tion of the quench, (b-c) Bin-averaged occupations of eigen- 
states for 20 realizations of the quench with bin size equal 
to N h = 20 ■ 32 = 640 in (b) and iV b = 20 • 128 = 2560 in 
(c). The bin averages are computed according to formulas 
(H2J and ■ The QMC prediction (|7J is represented by red 
solid line. The exponential thermal distribution with same 
average energy is represented by the dashed green line, (d) 
Inverse participation ratios of eigenstates given by Eq. (|14p for 
a single realization of H. 




FIG. 4: (Color online) Frames (a) and (b) have the same 
kind of content and the same plot legends as Figs. [3ja) and 
(b), respectively, but with the special diagonal element giving 
E av = —4. The large occupation of the ground state exhibited 
in frame (a) indicates the QMC condensation regime. The 
inset magnifies the small- (p) part of the main plot in (b). 
The density of states for this case is shown in Fig. [1] 



of random matrix. Overall, however, the averaged oc- 
cupation (pi)(E av ) shows quite a sharp transition to the 
condensation regime at an average energy of about -3.5. 
In all cases, the occupation of the lowest eigenstate lies 
significantly below the value given by Eq. ([§]) for a QMC 
ensemble in a macroscopic system with short-range in- 
teractions. 

In order to investigate the condensation regime further 
we introduce not one but three special diagonal elements 
with values -4, -4.5 and -5, and then choose the basis 
state corresponding to the diagonal element equal to -4 
as the initial state before the quench. The results are 
presented in Fig. [6] We observe that two eigenstates 
with the lowest and the second lowest energy are almost 
unoccupied. The third eigenstate from below has large 
occupation comparable to one, i. e. we observe conden- 
sation but not into the lowest energy state. The averaged 
occupations of the remaining eigenstates agree very ac- 
curately with the QMC formula. In the present case, 
choosing the basis states corresponding to the diagonal 
elements -5 or -4.5 as the initial states would lead to the 
condensation into the lowest or the second lowest eigen- 
states, respectively. We note here that condensation into 
the second or the third lowest energy state contradicts to 
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FIG. 5: (Color online) The onset of condensation into the 
ground state as a function of E av for the random matrix 
quench with one special diagonal element. Noisy blue line 
- occupation of the ground state pi for a single realization of 
random matrix with a given value of the special diagonal ele- 
ment. The solid red line shows p\ averaged over the random 
realizations. Dashed green line represents the prediction @ 
for the QMC ensemble in a macroscopic system. 




FIG. 6: Same as Fig.[4]but for the random matrix quench with 
three special diagonal elements -5, -4.5, -4 and E av = —4. 



the QMC statistics always predicting the condensation 
into the lowest energy state. We discuss this disagree- 
ment in subsection IIV Bl 
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FIG. 7: Same as Fig. SJb) but for the random matrix quench 
with non-zero entries in H selected according to the prescrip- 
tion for the Gaussian unitary ensemble as described in the 
text and with E av = —4. 



2. Sparse Gaussian unitary ensemble 

Here we investigate the case, when H is a random ma- 
trix of the same size and with the degree of sparseness 
as before but with the non-zero matrix elements selected 
according to the prescription for the Gaussian unitary en- 
semble: the real and the imaginary part of each non-zero 
off-diagonal element is chosen randomly from the Gaus- 
sian distribution with variance a = 1/ ^/2^ 1 the diagonal 
elements are real and chosen from the Gaussian distribu- 
tion with variance a = 1 . A diagonal element equal to -4 
is added, and the corresponding basis state is expanded 
in terms of cigenstates. The result presented in Fig. [7] 
exhibits a very good agreement with formula ([7]) for the 
QMC statistics. In the present case, the condensation 
regime was not reached for most random matrix realiza- 
tions, although a small minority of the realizations did 
exhibit it, implying the condensation precursor behavior 
similar to the one observed in Fig. [5] above E av = —3.5. 



B. Discussion: correlations, condensation and 
localization 

The numerical result presented in Fig. [2] indicates that 
for modest values of E av random matrix quenches lead to 
the ensembles that are approximately describable by the 
QMC statistics but with noticeable deviations from the 
prediction by Eq. 0. The QMC formulas © and 
correspond to the limit of N ^ 1 , which should indeed be 
appropriate for 4096 states involved. As discussed below, 
the deviations near the edges of the energy spectrum, can 
be attributed to the localization properties of the eigen- 
states involved. More difficult, however, is to understand 
the discrepancy in the middle range between E = — 4 and 
E = 4. Here we observe, that the deviation between the 
numerical result and formula ((7|) looks somewhat similar 
to the deviation between the same formula and the direct 
simulations of the QMC ensemble for the quantum sys- 
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terns with the number of levels of the order of 10 [24|,[25[. 
The results of those QMC simulations are not describ- 
able by formula (J7|) because of the finite- TV effects. We 
therefore suspect, that since the random matrix is sparse 
with each line having only about 30 non-zero off-diagonal 
elements, the basis states may form resonantly coupled 
groups. The number of these groups may indeed be of the 
order of 10. Within each group, the basis states dclocalize 
into a band-like structure, which is correlated in the sense 
that it is not obtained by a random rotation of the basis 
states belonging to the group. In the next order, these 
groups randomly couple to each other thereby forming 
the resulting eigenstates. If true, such a scenario sug- 
gests, that the eigenstates of the random-matrix Hamil- 
tonians considered emerge not through a truly random 
rotation of the entire 4096-dimensional Hilbcrt space, but 
rather through a rotation of a smaller subspace thereby 
leading to the deviations from the QMC predictions. The 
fact that band-like correlations lead to deviations from 
the QMC statistics will be further substantiated in Sec- 
tion [V] on disordered potentials. 

Figure [5] illustrates that random matrices with one spe- 
cial diagonal element exhibit a rather sharp transition to 
QMC-like condensation as a function of the special el- 
ement. In general, we observe in Figs. [3j SJ [6] and [3 
that the agreement with the QMC prediction improves 
as this prediction approaches the condensation regime. 
Once the condensation regime is reached, as in Figs. [4] 
and |6l the occupations of non-condensed states are very 
well describable by the QMC formula. 

The QMC ensemble implies condensation into the 
lowest energy state for macroscopic systems with finite 
Hilbcrt space, reasonably short range of interactions and 
almost any average energy per particle (in the spectral 
range of the system) . The density of states for the total 
energy of such a system has a narrow Gaussian maximum 
in the middle with half- width, which is much smaller than 
the entire range of the energy spectrum. In contrast, the 
density of states of random matrices considered in this 
work are approximately semicircle, i.e. the maximum in 
the middle is very broad. In this case, the QMC con- 
densation is still possible, but only for the values of E av 
sufficiently far from the middle of the spectrum. Another 
consequence of the absence of the narrow peak is the non- 
applicability of formula (|9]) for the condensed fraction 
(see Fig. [SJ. 

As shown in Fig. [6j when more than one special diag- 
onal element is introduced, the condensation occurs not 
necessarily into the ground state of the system, but pos- 
sibly into a higher energy eigenstate. This circumstance 
together with the fact that in Fig. [3jb) the deviations 
from the QMC statistics occur for the bin of the lowest- 
energy states indicates that the phenomenon of localiza- 
tion may be involved here. The QMC ensemble requires 
a fair sampling of the Hilbert space. This implies a broad 
participation of eigenstates ipi in a given basis state 4>k- 
But if some of the eigenstates are localized among very 
few basis states, they will not be present in most of the 



quench-induccd ensembles considered in this paper. This 
localization property can be characterized by the inverse 
participation ratio (IPR) of basis states (f>k in a given 
energy eigenstate ipi'- 

N 

IPR^£|(^ fc )| 4 (14) 
fc=l 

IPRi of the order of 1/N indicates that eigenstate ipi is 
delocalized; IPR 2 ^ 1/N indicates that the eigenstate is 
localized. We note here that the phenomenon of QMC 
condensation, if observed in the random matrix quench, 
indicates that the lowest-energy eigenstate is strongly lo- 
calized. 

In general, however, it is not just one eigenstate, but 
a group of eigenstates near both edges of the spectrum 
that are localized. A plot of the inverse participation 
ratio of eigenstates for a matrix ensemble with special 
element equal to -2 is presented in Fig. [3jd). The plot 
shows that states at very low energies are localized in the 
sense of having a rather high inverse participation ratio. 
Therefore, it is not surprising that the QMC statistics 
fails to describe their individual occupations. 

We define the "mobility edge" through the critical 
value IPR cr of the inverse participation ratio as 

IPR cr = 1.5 ■ IPRtyp, (15) 

where IPRt yp is the average value of the inverse participa- 
tion ratio of an eigenstate in the middle of the spectrum. 
We call an eigenstate localized, if its inverse participation 
ratio is larger than IPR cr . 

The idea is now to put all eigenstates below the mobil- 
ity edge into a single energy bin. If the QMC statistics 
works above the mobility edge, and all states below the 
mobility edge belong to a single bin, then the QMC pre- 
diction will be valid for the lowest energy bin as well. It 
is guaranteed by the normalization condition 

N 

(pi) = l-$>i>- (16) 

i=2 

The effectiveness of the above binning procedure can be 
observed by comparing Fig. G3h) and (c). In Fig. [3jb) 
the lowest energy bin is not large enough. Hence it does 
not include all localized eigenstates and exhibits clear 
deviation from the QMC prediction. In Fig. EHc), the 
bin size was increased, so that all eigenstates below the 
mobility edge belong to the first energy bin, which then 
exhibits a very good agreement with the QMC prediction. 

To summarize, our findings suggest that the QMC 
statistics works best, when the combined occupation of 
all eigenstates below the mobility edge is significant and 
all these states are grouped into a single energy bin. 

In the above discussion we neglected the existence of 
the second mobility edge near the top of the energy eigen- 
spectrum. In general, the occupations of eigenstates 
above the high-energy mobility edge also tend to devi- 
ate from the QMC statistics, but these occupations are 
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small when the condensation regime on the other end of 
the energy spectrum is approached. 



V. DISORDERED LATTICES 

In this section, we consider one particle in a single- 
band lattice. The Hamiltonian before the quench has 
disordered on-site energies £, and no hopping elements: 



(a) 



Ho = 



N 

i=i 



(17) 



Here is the particle annihilation operator for the zth 
site. The particle before the quench occupies the kth 
site and has energy E av — ej.. The quench consists of 
switching-on the nearest-neighbor hopping term 



n 



pert 



tij dj a j 



h.c. 



(18) 



where tij are the hopping elements, which are also dis- 
ordered. The notation implies the summation over 
nearest neighbors. 

The full Hamiltonian H = Ho + H pcr t is sparse, but, 
unlike the case of random matrices, it has non-zero off- 
diagonal elements located at positions that are not ran- 
dom. 

Physically, this situation may be realized with disor- 
dered potentials, when initially the potential barriers be- 
tween the wells are so large that the hopping can be ne- 
glected, and then the barriers are suddenly lowered, so 
that the particle can hop. Or, alternatively, a single par- 
ticle is injected into the fcth site of the lattice describablc 
by the full Hamiltonian H = Ho + H peI t ■ Both settings 
should be realizable in the experiments on optical lat- 
tices. They are also conceptually relevant to the experi- 
ments with disordered speckle potentials [lol liH 0j| 
and quasi-periodic potentials 1171 t hat reported the evi- 
dence of Anderson localization [2fJ ■ 
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A. Results and discussion 

1. Disordered lattices with real-valued hopping elements 

In this part, we investigate three lattices: 

(i) three dimensional (3D) face centered cubic (fee) lat- 
tice with 4000 sites; 

(ii) 3D simple cubic lattice with 4096 sites; and 

(hi) one-dimensional (ID) chain with 4096 lattice sites. 
For each lattice, we choose the on-site energies q ran- 
domly from a uniform distribution in the interval [-1,1] 
and hopping elements tij randomly from a uniform distri- 
bution in the interval [0,1]. The latter choice is motivated 
by the fact that for real potentials with deep enough 
wells, all hopping elements should be positive. The initial 



FIG. 8: (Color online) Quantum quenches for different disor- 
dered lattices with random on-site energies et selected from 
the interval [—1,1] and real random hopping elements ty se- 
lected from the interval [0, 1]. The type of lattice is indi- 
cated above each plot. The inset in (c) shows the results for 
the ID lattice with hopping elements selected from the in- 
terval [0.1, 1]. In all cases, the energy expectation value of 
the initial states is E av — —0.9. The numerical results rep- 
resent 600 realizations of each disordered lattice. To improve 
the statistics further, 20 basis states with energy expecta- 
tion value closest to -0.9 were included. The plots show the 
bin-averaged occupations of eigenstates (p) with bin sizes 
equal to 600-20-40 = 480000 in (a) and 600-20-32 = 384000 in 
(b) and (c). The solid red line represents the QMC prediction 
by Eq. (J7J). The dashed green line represents the exponential 
thermal distribution with same average energy. 
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site is characterized by the on-site energy E av w —0.9. 
The results are presented in Fig. [5] 

For all three lattices, both the QMC formula ([7]) and 
the exponential distribution can describe the occupa- 
tions of eigenstates only very roughly, with all numeri- 
cally generated distributions exhibiting large deviations. 
These deviations may be related to band-like correlations 
mentioned in subsection IIVB1 because the densities of 
states for all three lattices considered have traces of the 
band structures associated with the limit of no onsitc 
disorder and translationally invariant values of t irj . In 
the translationally-invariant case, the occupation of each 
eigenstate (Bloch state) would have equal value, which is 
clearly different from the QMC prediction. On the other 
hand, we note, that the deviations from the QMC statis- 
tics do not appear to have systematic character, which 
implies that the approximation based on the QMC en- 
semble is, probably, the best one can do without perform- 
ing a system-specific analysis. In particular, the QMC 
statistics can be reasonably used to estimate the proba- 
bility of a particle being found above or below the mo- 
bility edge after the particle is injected into a site with a 
given value of the on-site energy E av . 

A priori, we expected that the result for the fee lattice 
would be closest to the random matrix case because of 
the largest number of nearest neighbors (equal to 12) and 
thus closest to the QMC prediction. As seen in Fig. [SJa) , 
this expectation was indeed confirmed. 

On the opposite end is the result shown in Fig. [HJc) 
for ID lattice, which exhibits sharp edge and a peak, 
both close to E av . The above peak can be understood by 
observing that each site of the ID lattice has only two 
nearest neighbors. The corresponding hopping elements 
are chosen randomly from the interval [0,1], which means 
that there is a non-zero chance for both of them to be 
very close to zero. In such a case, the initial state is close 
to a strongly localized eigenstate, which means that the 
occupation of a single eigenstate is comparable to 1. The 
inset of Fig. |5[c) shows the result for a simulation with 
hopping elements chosen from the interval [0.1,1], i.e. the 
possibility of the near-zero hopping elements is excluded. 
The peak at the average energy vanishes, as expected 
from the above qualitative picture. For the other two 
lattices in Fig. the probability of all hopping elements 
around a given site to be close to zero is suppressed due 
to higher numbers of nearest neighbors. Hence the peak 
disappears. 



2. Disordered lattices with complex-valued hopping elements 

In experiments with ultra cold atoms on optical lat- 
tices, it may be possible to break the time-reversibility 
by simulating the effect of a staggered magnetic field act- 
ing on a charged particle [26|. In this case, the hopping 
elements tu are complex. We investigated this case by 
studying the fee lattice with complex-valued hopping el- 
ements. 




0.0000 1 



FIG. 9: (Color online) Same as Fig.[8ja) but with the on-site 
energies and the hopping elements selected according to the 
prescription for the Gaussian unitary ensemble as described 
in the text. Energy expectation value of the initial state is 
E av = —3. 600 independent realizations of the disordered 
fee lattice were used. The size of the energy bin is iVj, = 
600 ■ 40 = 24000. The inset magnifies the small- (p) part of 
the main frame. 



In this case, we chose the on-site energies from a Gaus- 
sian distribution with variance a = 1, and both Re(ijj) 
and lm(tij) (independently) from a Gaussian distribu- 
tion with variance a = 1/ \/2. Initially, the particle is 
located on a site with E av = —3. The result is presented 
in Fig. IH1 In contrast to the results for the lattices with 
real and positive hopping elements, the agreement of the 
bin-averaged occupations of eigenstates with the QMC 
formula ([7]) is very good. This is consistent with the ex- 
pectation that, in the present case, band-like correlations 
are destroyed more efficiently due to the random phases 
of the hopping elements. Also E av in this case is closer 
to the QMC condensation threshold. 



VI. CONCLUSIONS 

We have investigated quantum ensembles emerging as 
a result of off-diagonal quenches of sparse random-matrix 
Hamiltonians, and disordered potentials. The resulting 
ensembles arc controlled by the energy expectation value 
of the initial quantum state. We have introduced a new 
class of random-matrix Hamiltonians with one special di- 
agonal element and shown that quantum quenches with 
these Hamiltonians exhibit transition to the QMC con- 
densation. In general, random-matrix eigenstates close 
the edges of the energy spectra tend to be more local- 
ized in the original basis, and hence the occupations of 
these eigenstates follow the QMC statistics not so well. 
However, if a mobility edge for the eigenstates is intro- 
duced, and all states below the mobility edge are grouped 
into a single energy bin, then the QMC statistics is recov- 
ered. In the case of disordered potentials, our simulations 
indicate that the resulting ensembles, in general, do not 
exhibit quantitative universality. Nevertheless, the QMC 
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statistics can still serve at least as a crude approximation. when the number of the nearest neighbors on the lattice 
The quality of this approximation improves significantly is increased, or when the time-reversibility is broken. 
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